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Summary 

This  report  along  with  a  video  demonstration  tape  documents  the 
efforts  expended  In  Investigations  conducted  at  Magnavox  Into  the 
development  of  real-time  algorithm  techniques  which  can  be  used  In 
FLIR  Imagery  for  the  detection  of  cables  In  the  field  of  view. 

A  combination  of  "seml-llnear"  local  line  detectors  and  ’’semi- 
local”  line  discriminators  were  developed  to  optimize  detection  and 
false  alarm  rejection. 

Conclusions  and  recommendations  for  future  efforts  ^n  this  area 
are  discussed  In  Section  6. 
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1.0 


INTRODUCTION 


This  final  technical  report  covers  the  progress  and  efforts  expended 
by  Magnavox  Government  and  Industrial  Electronics  Company  (MAGIEC)  towards  the 
objectives  of  Contract  DAAK70-81-C-0177  for  the  period  of  1  October  1981  to 
1  September  1982,  and  Is  submitted  together  with  a  video  demonstration  tape  In 
compliance  with  Contract  Data  requirements  List  Item  AOOl.  This  contract 
comprises  an  Investigation  Into  the  real-time  algorithm  techniques  which  may  be 
used  In  FLIR  Imagery  for  the  detection  of  cables  In  the  field  of  view.  The 
video  tape  shows  examples  of  the  developed  algorithm  applied  to  video  scenes 
containing  cables. 

1.1  BACKGROUND 

It  has  long  been  a  problem  for  pilots  of  low  flying  helicopters  to 
detect  and  avoid  wire  obstacles  such  as  telephone  lines  and  power  cables.  Past 
experience  has  shown  that  wire  strikes  may  occur  during  routine  day  operation. 
With  Increased  emphasis  on  Nap-of-the-Earth  (NOE)  flight  as  a  tactic  to  avoid 
detection  by  threat  radars  and  optical  systems,  the  wire  avoidance  problem 
becomes  even  more  critical.  In  addition,  NOE  operations  are  now  expanding  Into 
periods  of  limited  visibility  and  darkness  through  the  use  of  IR  systems  such  as 
the  Pilot's  Night  Vision  System  (PNVS)  on  the  Advanced  Attack  Helicopter  (AAH). 

NOE  flight  under  the  best  of  conditions  Is  a  demanding  pilot  task. 

NOE  flight  at  night  with  a  PNVS  creates  an  environment  In  which  the  pilot  must 
work  at  near  maximum  capability.  It  Is  In  this  environment  that  the  pilot  Is 
easily  detracted  by  other  Imagery  In  the  scene  or  other  flying  tasks  and  the 
ability  to  perceive  faint  wire  obstacles  Is  seriously  degraded.  Therefore,  the 
problem  addressed  In  this  contract  Is  to  determine  methods  which  will  detect 
wire  obstacles  In  the  scene  video  of  a  PNVS. 

Previous  solutions  to  the  wire  avoidance  problem  have  addressed  the 


general  case  of  both  day/nlght  operations  and  have  been  Independent  of  other 
on-board  sensors.  These  solutions  have  required  the  addition  of  dedicated 
sensors  to  the  helicopter  resulting  in  greatly  increased  weight,  complexity,  and 
cost.  In  addition,  the  effectiveness  of  these  sensors  has  been  marginal.  For 
the  AAH,  however,  the  PNVS  represents  a  high  performance  FLIR  which  through  more 
effective  utilization  can  provide  imaging  of  wire  obstacles.  The  solution 
investigated  here  is  the  extraction  of  wire  obstacles  from  video  Imagery  using 
advanced  image  processing  techniques.  With  the  development  of  the  common  module 
Digital  Scan  Converter  (DSC),  most  of  the  required  hardware  will  already  be 
on-board  the  helicopter.  This  approach  also  has  the  advantage  of  being  passive 
whereas  the  previous  solutions  have  required  active  systems,  thus  Increasing  the 
helicopter's  vulnerability  In  a  combat  environment. 


2.0 


THE  PROBLEM 


A  capability  for  cable  avoidance  through  passive  IR  image  processing 
techniques  requires  the  ability  to  discern  thin  lines  in  the  presence  of 
competing  scene  structure,  fixed  or  moving  pattern  noise,  and  random  noise, 
based  solely  on  the  brightness  distribution  measured  over  the  scene.  This  is  a 
difficult  requirement  under  the  best  of  circumstances,  and  necessitates  an 
algorithm  sensitive  to  those  characteristics  in  the  brightness  distribution 
which  are  unique  to  lines. 

Such  an  algorithm  must  be  able  to  distinguish  between  lines  and  scene 
clutter  such  as  edges,  points,  clusters,  and  similar  features  (intensity 
variations)  which  arise  from  buildings,  roads,  foliage,  etc.,  found  in  any 
terrain  scene.  These  features  can  effectively  hide  or  break  up  the  unobtrusive 
signatures  of  hanging  wires. 

An  even  more  diffucult  task  is  for  the  algorithm  to  distinguish 
between  lines  resulting  from  hanging  cables  and  lines  due  to  other  phenomena. 
The  most  serious  problems  arise  from  streaking  caused  by  detector 
non-uniformities  and  from  raster  crawl  caused  by  Incompatibilities  in  FLIR  and 
TV  scanning  modes.  However,  the  advent  of  the  digital  scan  converter 
technology,  with  its  TV  compatible  output  and  its  gain  and  level  equalization 
preprocessing  will  considerably  reduce  these  problems.  The  algorithm  should  be 
able  to  discriminate  against  any  residual  effects.  Lines  due  to  other  scene 
components  besides  cables  may  not  be  so  easily  dealt  with,  but  on  the  other 
hand,  are  not  that  severe  a  problem.  Nevertheless,  the  ideal  algorithm  should 
be  able  to  distinguish  between  these  lines. 

Random  system  noise  and  high  spatial  frequency  scene  fluctuations  due 
to  foliage,  etc.,  has  small  line  segment  structure  and  point  structure.  The 
line  segment  structure  results  in  algorithm  false  alarms,  while  point  structure 


near  true  lines  alters  the  line  signatures  and  results  In  algorithm  failure  to 
detect  real  signals.  Digitization  noise,  and  signal  aliasing  effects  of  the 
sampling  process  Intrinsic  to  digitizing  a  scene  also  have  a  corrupting  effect 
on  line  signatures.  A  good  algorithm  must  be  able  to  operate  effectively  In 
poor  slgnal-to-nolse  ratio  environments. 

The  problem  of  cable  avoidance  from  an  algorithm  point  of  view  can  be 
separated  Into  two  distinct  tasks:  (1)  detecting  potential  line  (cable)  segments 
In  the  scene  In  the  presence  of  noise  and  clutter  and  (2)  connecting  the  proper 
detected  line  segments  Into  extended  lines  further  eliminating  false  alarms. 

The  first  task  requires  a  local  algorithm  utilizing  an  operator  that  Is 
sensitive  to  the  scene  Intensity  statistics  In  the  neighborhood  of  a  given  scene 
element,  or  pixel*.  Its  difficulty  Is  In  discriminating  between  a  line  element 
and  a  non-line  element  (e.g.,  an  edge  element  or  point  element).  An  overly 
restrictive  algorithm  will  discard  too  many  real  line  segments  and  an  overly 
loose  algorithm  will  keep  too  many  clutter  elements  (false  alarms).  The  second 
task  requires  a  global  algorithm  with  the  capability  for  remembering  all 
detected  line  segments  from  across  the  entire  scene,  discriminating  between 
potential  cable  segments  and  false  alarms,  and  connecting  the  cable  segments 
Into  extended  lines  which  are  true  representations  of  the  original  cables.  This 
algorithm  may  require  bridging  large  gaps  caused  by  Intervening  scene  objects. 

There  are  two  constraints  that  make  the  problem  even  more  difficult. 
First,  the  complete  algorithm  must  operate  In  nearly  real  time,  l.e..  In  a  time 
short  enough  to  follow  scene  movement  Incurred  due  to  vehicle  or  FLIR  motion, 
and  In  any  case  short  enough  to  allow  the  pilot  to  take  evasive  action.  Second, 
the  algorithm  must  be  Implementable  In  a  small,  light  weight  format 

*A  picture  element  Is  completely  determined  by  Its  position  coordinates  and  Its 
brightness  value. 


(microprocessor,  etc.)  for  Inclusion  In  the  on-board  electronics. 

Due  to  the  limited  scope  of  this  effort  It  was  decided  to  confine  the 
Investigation  to  the  local  detection  algorithm  since  this  appears  to  be  the 
critical  phase  of  the  problem.  It  may  be  possible  to  allow  the  eye  to  connect 
line  segments  that  have  been  sufficiently  enhanced  to  be  Immediately  obvious  to 
the  pilot.  Furthermore,  failures  of  the  detection  algorithm  obviates  the  need 
for  a  connection  algorithm.  However,  some  general  concepts  and  aspects  of  a 
connection  algorithm  will  be  discussed  In  terms  of  the  future  directions  the 
algorithm  should  take  for  a  successful  completion  of  the  task. 


I 


THE  CANDIDATE  ALGORITHMS  FOR  LOCAL  DETECTION 


Any  algorithm  that  detects  lines  must  be  based  on  some  essential 
characteristic  or  group  of  characteristics  of  a  line  that  make  It  different  from 
other  picture  elements.  The  essential  Ingredient  of  a  linear  feature  Is  that  It 
have  continuity  In  one  direction,  but  represent  some  deviation  from  the  norm  In 
any  other  direction.  For  lines.  It  Is  the  picture  Intensity  that  Is  continuous 
and  slowly  varying  along  the  line  (l.e.,  on  and  In  the  direction  of  the  line), 
but  has  an  extremum  In  any  other  direction.  Thus,  a  line  has  three  key 
qualities:  , 

1)  A  direction 

2)  A  finite  extent  In  that  direction 

3)  An  extreme  Intensity  value  across  that  direction.  (For  real  lines 
of  finite  width  the  Intensity  gradient  should  be  a  maximum  along  the 
perpendicular  to  the  line's  direction.)  The  candidate  detection  algorithms  will 
make  use  of  some  or  all  of  these  qualities. 

3.1  MEDIAN  FILTERING 

One  algorithm  which  makes  use  of  a  non-unique  statistical  feature  of  a 
line  applies  the  compliment  of  a  median  filter  to  the  scene.  Since  the  elements 
of  a  line  are  greater  (or  less)  than  their  off  line  neighbors,  they  will  never 
be  the  median  value  In  their  neighborhood,  and  In  fact  these  line  elements  make 
a  maximum  difference  with  the  neighborhood  median.  Consequently,  If  we  subtract 
the  neighborhood  median  from  each  pixel,  we  should  get  extrema  on  lines  and  a 
washed  out  scene  away  from  lines.  However,  while  this  technique  can 
discriminate  between  lines  and  edges.  It  cannot  do  so  between  lines  and  points 
since  point  elements  by  definition  are  greater  (or  less)  than  their  neighbors. 
Because  this  algorithm  falls  to  take  Into  account  the  directivity  and 
connectedness  of  lines  vis  a  vis  points.  It  enhances  points  as  well  as  lines. 


The  resulting  ‘'line-enhanced"  scene  is  much  noisier  than  the  original  scene. 
Although  originally  it  was  planned  to  pursue  this  technique  based  on  our 
previous  observations  of  its  line  enhancing  capability,  it  has  been  discarded 
because  of  its  noise  enhancing  side  effect.  Instead  a  search  was  undertaken  for 
a  technique  which  makes  use  of  three  key  qualities  of  a  line  mentic'ed  above. 

3.2  LOCAL  DETECTOR  ALGORITHM  STRUCTURE 

A  line  element  has  a  unique  relationship  to  its  neighbors  (l.e.,  to 
other  nearby  line  elements  and  to  neably  off-line  elements)  determined  by  those 
qualities  which  define  a  line.  To  detect  a  line  element  we  must  first  determine 
the  value  of  the  neighboring  elements  or  at  least  sample  the  neighborhood,  and 
then  verify  that  the  proper  relationship  exists  for  some  choice  of  direction 
(l.e.,  for  some  assignment  of  on-line  and  off-line  elements).  This  is 
accomplished  by  choosing  a  set  of  sampling  templates,  each  member  of  which 
assumes  a  different  orientation  of  the  line  and  assigns  the  identity  of  the 
on-line  samples  and  the  off-line  samples  which  will  be  tested,  based  on  the 
assumed  direction.  If  each  member  sampling  template  is  applied  to  a  picture 
element  and  its  neighborhood,  the  relationship  of  the  picture  element  to  the 
neighborhood  samples  will  conform  to  the  line  detection  criteria  only  when  the 
picture  element  is  Indeed  a  line  element  and  then  only  when  the  assumed 
direction  is  correct.  The  algorithm  can  then  be  structured  as  follows: 

a.  Apply  the  sampling  template  set  to  each  element  and  its 
neighborhood. 

b.  Apply  detection  criteria  for  evidence  of  a  line  for  each  template. 

c.  Compare  evidence  for  each  direction  and  choose  direction  which 
gives  "best”  evidence  for  a  line,  if  one  exists. 

d.  If  no  evidence  exists,  or  evidence  is  not  sufficient  (l.e.,  below 
threshold)  assume  no  detection. 


3.2.1 


The  Sampling  Template  Set 


The  size  and  sampling  configuration  of  the  sampling  template  set  are 
Important  In  determining  the  direction  and  width  of  a  line.  A  3  x  3  sampling 
template  set  allows  only  four  (4)  different  directions  giving  a  resolution  of 
45”.  A  5  X  5  set  allows  for  eight  (8)  possible  directions  (a  resolution  of 
22.5”)  although  the  sampling  configuration  Is  somewhat  awkward  for  some 
directions.  Figure  3.1  Is  an  example  of  a  3  x  3  neighborhood  sampling  template 
for  which  the  assumed  line  direction  Is  vertical.  82  Is  the  queried  pixel,  Bi 
and  83  are  Its  neighboring  on-line  pixels,  and  the  A's  and  C's  are  the  parallel 
off-line  neighbors. 
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Figure  3.1 


80th  the  3x3  and  5x5  sets  require  that  the  samples  be  contiguous.  This 
means  that  only  thin  lines  (approximately  1  pixel  wide)  will  be  detected.  If 
the  samples  were  separated,  thicker  lines  would  also  be  detected.  In  order  to 
have  a  template  set  with  at  least  8  possible  orientations  (22.5”  resolution)  and 
with  a  non-contlguous  symmetric  sampling  configuration,  a  9  x  9  neighborhood 
need  be  sampled.  Therefore  a  9  x  9  sampling  template  set  was  chosen  for  this 
Investigation  and  Is  shown  In  Figure  3.2. 


3.2.2  LOCAL  LINE  DETECTION  CRITERIA 

To  constitute  a  line  detection  the  on-line  sample  elements  must  be 
greater  (or  less)  than  the  neighboring  off-line  sample  elements  to  either  side. 
Mathematically,  this  can  be  formulated  by  computing  the  differences  of  the 
on-line  sample  elements  minus  the  off-line  sample  elements  (e.g.,  Bi-A^,  Bi-Ci) 
and  requiring  that  they  all  have  the  same  sign.  However,  there  are  various  ways 
to  formulate  this  criteria,  each  with  a  somewhat  differing  result.  These 
approaches  can  be  classified  as  linear,  seml-llnear,  and  non-linear.  They 
differ  only  in  the  degree  to  which  they  employ  the  unique  qualities  of  a  line. 

3 . 2 . 2 . 1  The  Linear  Algorithm 

The  linear  approach  requires  that  the  average  of  the  on-line  samples 
be  greater  (or  less)  than  the  average  of  all  the  off-line  samples  1^  some 
threshold  amount  or  no  detection  is  encountered.  If  a  ■  £Ai,  b  "  ZBi,  and 
c  ■  rCi  (ref.  Figure  3.2)  and  if  E  is  the  evidence  for  a  line  then  the  linear 
algorithm  can  be  simply  stated  mathematicaly  by  the  following  equations, 

E  -  2b  -  a  -  c  ,  if  Ie]  >  t;  (3.1) 

E  >  0  ,  otherwise. 

Due  to  the  non-uniqueness  of  this  statistical  property,  this  method  responds  to 
edges  and  Isolated  points  as  well  as  to  lines  although  its  line  response  is  much 
stronger.  On  the  other  hand,  because  it  averages  over  all  off-line  elements  at 
once,  it  tends  to  filter  out  the  corruptive  effect  of  noise  on  line  segments. 
Nevertheless  it  was  decided  to  eliminate  the  linear  algorithm  Immediately,  based 
on  its  inability  to  discriminate  against  edges  and  points. 

3. 2. 2. 2  The  Non-Linear  Algorithm 


2 


if  Bi  -  Ai  t  and  -  C£  >  t  for  all  1  (3.2) 

or  If  Bf  -  Ai  £  t  and  Bj  -  Cj  £  t  for  all  1; 

E  "  0,  otherwise 

This  algorithm  discriminates  against  both  edges  and  Isolated  points,  but  it  has 
a  problem  dealing  with  short  gaps  in  a  line  caused  Interfering  features  and 
(especially)  corruptive  noise.  Vhlle  this  algorithm  is  effective  in  eliminating 
false  alarms,  it  has  a  low  detection  capability  for  lines  in  the  presence  of 
noise.  The  above  algorithm  was  coded  and  tested  along  with  the  seml-llnear 
algorithm. 

3. 2. 2. 3  The  Seml-Llnear  Algorithm 


The  semi-linear  algorithm  is  a  compromise  between  the  linear  and  the 
non-linear  algorithms.  While  it  does  not  average  over  all  off-line  samples, 
neither  does  it  treat  each  off-line  sample  separately.  Rather,  it  averages 
separately  over  the  two  sets  of  off-line  samples  (one  set  to  either  side  of  the 
assumed  line).  That  is,  the  seml-llnear  algorithm  requires  that  the  average  of 
the  on-line  samples  be  greater  (or  less)  than  either  off-line  set-average  ^ 
some  threshold  amount.  An  example  of  a  seml-llnear  detection  algorithm  is  given 


by  the  equations. 


E  -  2b  -  a  -  c 


if  b  -  a  >  t  and  b  -  c  >  t 


or  if  b  -  a  <  t  and  b  -  c  <  t; 


(3.3) 


E  ■  0  ,  otherwise. 

This  seml-llnear  approach  eliminates  edges  but  is  still  somewhat 


responsive  to  isolated  points,  although  less  so  than  the  linear  algorithm. 
However,  it  is  less  Influenced  by  corrupting  noise  near  line  elements.  This 
algorithm  was  also  coded  and  evaluated  in  this  study. 


3. 2. 2. 4  A  Variation 


One  can  choose  another  estlaate  for  the  evidence  besides  the  sum  of 
differences  given  In  equations  3«1>  3.2,  and  3.3.  For  example,  In  the  linear 
approach  one  could  choose  the  product  of  the  difference, 

E  -  (b-a)  •  (b-c)  If  (b-a)  >  t  and  (b-c)  >  t 

or  If  (b-a)  <  t  and  (b-c)  <  t;  (3.4) 

E  -  0  ,  otherwise. 

This  variation  should  Improve  the  discrimination  between  lines  and  edges.  This 
Is  because  the  product  favors  lines  having  similar  off-line  statistics  to  either 
side  of  the  line,  vis  a  vis,  lines  whose  off-line  samples  to  one  side  of  the 
line  are  much  stronger  than  the  off-line  samples  to  the  other  side  of  the  line, 
a  condition  which  strongly  resembles  an  edge  corrupted  by  noise. 

3.2.3  MAGNITUDE  AND  DIRECTION  ALGORITHMS 

Once  the  evidence  for  a  line  has  been  computed  for  each  orientation  of 
the  sampling  set  It  then  needs  to  be  evaluated  to  determine  which  orientation 
has  the  strongest  evidence.  Since  the  gradient  Is  greatest  when  traversing 
across  the  line  perpendicularly,  the  evidence  (differences)  should  be  strongest 
when  the  chosen  on-line  elements  are  most  closely  aligned  with  the  actual 
direction  of  the  line.  If  we  take  as  the  direction  of  the  line,  the  direction 
associated  with  the  sampling  template  which  has  the  strongest  evidence,  we 
should  be  able  to  evaluate  the  direction  of  the  actual  line  to  within  +  11.25”. 
At  this  point  we  can  also  define  some  measure  of  the  strength  of  the  line  based 
on  the  evidence. 

3. 2. 3.1  The  Maximum  Evidence  Approach 

The  most  straightforward  way  Is  to  Just  choose  the  value  of  the 
maximum  evidence.  If  VAL  Is  the  magnitude  of  the  line,  and  THETA  its  direction 
(l.e.,  the  angle  It  makes  with  the  horizontal  axis  as  shown  In  Figure  3.3),  then 


VAL  -  max  {E^J  -  Efca,  (3.5 
where  Ei^  Is  the  evidence  from  the  kth  sampling  template  and  km  is  the  index  of 
the  template  having  the  maximum  evidence.  The  angle  associated  with  the  k^^ 


template  is  given  by. 


(k-1)  IT  .. 

■  — g  ■ '  radians 


Figure  3.3 

VAL  can  be  thresholded  to  reduce  false  alarms.  This  algorithm  was  tested  and 
evaluated  In  the  study. 

3. 2. 3. 2  The  Statistical  Approach 

Another  more  complex  algorithm  uses  circular  statistics  to  develop  a 
more  accurate  estimate  of  the  direction  and  strength  of  a  line.  The  direction 
can  be  obtained  from  the  weighted  mean  orientation  0  by 

X  -  Z  Ek  cos  (2  01c)  (3.5 

Y  -  Z  Ek  sin  (2  ©ic)  (3.6 

0  -  (1/2)  tan  -1  (Y/X).  (3.7 

THETA  >  ©Ic  nearest  to  9.  (3.8 

The  strength  of  the  line  can  be  derived  from  the  variance  of  the  weighted 
distribution  of  directions  from  the  mean  direction.  If  the  variance  is  given 


p2  -  1  -  (X2  +  y2)/(i;  Eic)2 
then  VAL  -  1  -  p2  -  (X2  +  y2)/z  Eic)2 


VAL  Is  now  a  measure  of  the  unanimity  or  agreement  of  the  statistical 
directional  data  with  the  weighted  mean  orientation* 

3. 2. 3. 3  The  Symmetry  Parameter 

In  addition  one  can  define  a  third  parameter  which  can  be  used  as  a 
figure  of  merit  of  a  line.  This  parameter  depends  on  the  symmetry  of  the  line 
across  the  direction  of  the  line,  and  Is  given  by. 


MU  -  (2b  -  a  +  c)/4 

(3.11) 

(SIGMA)2  - 

[2(b  -  MU)2  +  (a  -  MU)2  +  (C  -  MU)2]/4 

(3.12) 

KAPPA  -  (2b  -  a  -  c)/4 

(3.13) 

and 

_  KAPPA 

SIGMA 

(3.14) 

where 

0  <  SYM  <  1. 

SYM  Is  a  measure  of  the  degree  to  which  the  off-'llne  samples  to  either  side  of 
the  line  have  the  same  value.  As  such  this  parameter  helps  discriminate  against 
edge  like  features.  Often  times  an  edge  which  is  corrupted  by  noise  can  meet 
the  statistical  requirements  of  a  line,  although  Its  symmetry  (as  measured  by 
SYM)  would  be  low.  SYM  can  also  be  thresholded  to  reduce  false  alarms.  An 
algorithm  using  THETA,  VAL,  and  SYM  to  evaluate  the  evidence  for  a  line  was  used 


In  this  study 


4.0 


THE  SEMI-LOCAL  LINE  DISCRIMINATOR 


In  order  to  reduce  the  false  alara  rate  without  overly  degrading  the 
detection  capability,  a  discrimination  technique  Is  required  that  Is  Independent 
of  the  detection  mechanism.  The  previous  algorithms  made  use  of  the  local 
property  of  a  line  that  the  pixel  value  Is  an  extremum  on  a  line  when  traversing 
across  the  line  and  that  the  extremum  Is  largest  when  the  direction  traversed  Is 
normal  to  the  line  direction  (Property  3  of  section  3.0).  What  has  not  been 
employed  Is  Property  2  of  section  3.0,  namely,  that  a  line  has  finite  extent 
along  Its  direction.  That  is,  there  Is  some  continuity  along  the  direction  of 
the  line  In  both  magnitude  and  direction,  particularly  In  direction  since 
fluctuations  In  background  (off-axis  texture,  etc.)  could  affect  the 
magnitude.  Therefore,  a  detected  line  segment  can  be  expected  to  have  on-line 
neighbors  which  are  also  detected  and  which  have  the  same  direction.  To  make 
this  determination  about  a  detected  line  segment  requires  extending  the 
neighborhood  of  consideration,  l.e.,  merging  neighborhoods  In  the  direction  of 
the  detected  line  segment.  In  this  sense  It  Is  a  semi-local  algorithm,  a 
compromise  between  the  purely  local  detection  algorithm  and  the  global 
connection  algorithm  which  would  link  locally  deleted  line  sections  Into 
continuous  lines  which  extend  across  the  scene. 

The  semi-local  line  discriminator,  examines  a  detected  line  element  to 
see  If  has  extent  along  Its  direction,  accepting  only  those  detected  line 
elements  which  have  detected  line  elements  at  either  end  of  Its  direction,  or  at 
least  Its  Immediate  neighborhood  which  line  segments  have  the  same  or  similar 
direction  and  perhaps  the  same  sign.  (The  sign  of  a  line  element  will  be 
considered  positive  If  the  local  on-line  extremum  Is  a  maximum,  and  negative  If 
It  Is  a  minimum.  Having  the  same  sign  Is  a  compromise  between  having  the  same 
magnitude  and  requiring  no  magnitude  correlation  whatever. )  To  make  these 


comparisons  requires  more  memory  since  It  becomes  necessary  to  store  the 
direction  and  magnitude  (or  just  sign)  Information  of  the  neighboring  detected 
line  elements  for  future  comparisons.  There  are  numerous  variations  possible  In 
this  correlation  process.  Two  algorithms  were  chosen  for  this  study,  one  based 
on  nearest  neighbors  and  one  based  on  second  nearest  neighbors. 

4 . 1  NEAREST-NEIGHBOR  CORRELATION 

The  first  algorithm  that  was  used  tested  only  nearest  neighbors  for  a 
correlation  of  direction  and  sign.  For  nearest  neighbors  there  are  only  four 
orientations  where  a  direct  extension  of  the  detected  line  element  leads  to  an 
unambiguous  choice  of  neighboring  on-line  elements  (and  their  orientation). 

Since  there  are  eight  possible  line  orientations,  the  algorithm  used  did  not 
require  specific  neighboring  elements  to  be  detected  line  elements  of  a  given 
orientation  and  sign.  Instead  It  was  required  that  the  neighborhood  contain  at 
least  a  specified  number  (N)  of  detected  line  elements  of  the  same  orientation 
and  sign.  To  carry  out  the  comparison,  each  detected  line  element  was  assigned 
a  direction  number  between  1  and  8  and  a  sign,  positive  for  a  bright  line  and 
negative  for  a  dark  line.  (For  example,  a  bright  vertical  line  Is  assigned  a 
value  of  +5  whereas  a  dark  line  at  45”  Is  assigned  a  value  of  -3.)  If  N  Is  set 
to  2  then  Figure  4.1  shows  sample  3X3  neighborhoods  of  detected  line  elements 
for  which  the  central  element  Is  considered  a  valid  line  element. 

4.2  SECOND-NEAREST  NEIGHBOR  CORRELATION 

If  second  nearest  neighbors  are  used,  the  test  Is  more  stringent  since 
It  requires  a  larger  extent  of  the  continuity  of  the  line  segment.  Furthermore, 
all  eight  possible  line  orientations  lead  to  an  unambiguous  choice  of  which 
second  nearest  neighbors  should  also  be  line  segments  and  what  their 
orientations  should  be.  An  extrapolation  of  the  central  detected  line  element 
In  a  5  X  5  neighborhood  should  Intersect  two  second  nearest  neighbors.  The 
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algorithm  used  in  this  study  then  requires  that  the  two  Intersected 
second-nearest  neighbors  also  be  detected  line  elements  with  the  same 
orientation  and  sign  as  the  central  detected  line  element.  No  requirements  are 
specified  for  the  other  nearest  neighbors.  Figure  4.2  shows  sample  second 
nearest  neighbor  configurations  for  which  Che  central  element  Is  considered  a 
valid  line  element. 

4 . 3  EVALUATION  OF  THE  SEMI-LOCAL  LINE  DISCRIMINATOR 

The  algorithm  using  nearest  neighbors  was  coded  and  tested  In  such  a 
way  that  one  could  vary  N,  the  minimum  number  of  nearest  neighbors  required  to 
be  line  elements  of  the  same  orientation  and  sign.  This  algorithm  turned  out  to 
be  quite  good,  yielding  a  significant  Improvement  In  the  false  alarm 
rejection.  The  best  case  was  for  N  equal  to  2,  and  this  value  was  used  In  all 
subsequent  test.  However,  there  was  still  room  for  more  Improvement. 

The  algorithm  using  second  nearest  neighbors  was  also  coded  and 
tested.  The  result  was  a  remarkable  Improvement  In  the  rejection  of  false 
alarms,  even  better  than  the  algorithm  using  only  nearest  neighbors. 

Since  the  above  two  semi-local  line  discriminators  were  Independent, 
they  could  effectively  be  Joined  to  yield  a  combined  algorithm  that  was  better 
than  either  alone.  This  was  tried  and  found  to  be  the  case.  In  fact  the 
Improvement  over  the  local  detector  results  were  so  dramatic  that  any 
Improvements  resulting  from  Interchanging  the  various  local  detector  algorithms 
was  comparatively  unimportant.  The  combined  semi-local  discrimination  algorithm 
was  tried  with  each  of  the  potential  local  line  detection  algorithms  and  found 
to  work  approximately  equally  well  with  each.  The  key  consideration  In  choosing 
the  local  detection  algorithm  Is  therefore  the  computation  time.  Choosing  the 
fastest  local  algorithm,  namely  one  which  takes  the  evidence  as  a  sum  of  on-line 
minus  off-line  differences  Instead  of  a  product,  and  one  which  simply  chooses 
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Figure  4.2  Sample  Second  Nearest  Neighbor  Configurations  Yielding  a 
Valid  Central  Line  Element 


the  ■axlnum  evidence  as  Indicating  the  line  direction,  should  yield  the  best 
overall  result.  Since  there  Is  no  significant  difference  In  computation  time 
between  the  non-linear  and  the  seml-llnear  local  detection  algorithms,  the 
preferred  seml-llnear  algorithm  was  chosen. 

The  possibility  of  combining  the  local  detection  and  semi-local 
discrimination  algorithms  with  the  median  filtering  algorithms  was  considered. 
This  was  tested  by  first  subtracting  the  median  filtered  scene  from  the  original 
scene  to  give  the  "enhanced"  scene.  Then  the  algorithms  evaluated  above  were 
applied  to  this  scene.  The  results  were  appreciably  degraded  from  the  results 
without  the  median  filter  algorithm,  due  principally  to  the  Increased  false 
alarms  caused  by  the  point  enhancing  property  of  the  median  filter  algorithm. 

One  result  of  the  Improved  discrimination  of  the  semi-local  algorithm 
Is  that  the  detection  thresholds  can  be  lowered  substantially  without  a  large 
Increase  In  false  alarms.  The  reduction  In  optimum  threshold  yields  a  more 
sensitive  total  algorithm.  However,  Its  sensitivity  is  still  somewhat  limited. 
What  Improvements  can  be  made  In  sensitivity  In  subsequent  efforts  Is  discussed 


5.0  FORMAT  OF  THE  SOFTWARE  IMPLEMENTATION 

In  order  to  establish  a  framework  in  which  to  compare  the  variations 
In  the  cable  enhancement  algorithm.  It  was  decided  to  write  a  generalized  main 
routine  which  would  access  the  various  parts  of  the  algorithm  and  their 
variations  on  a  subroutine  call  basis.  In  this  way  the  different  sub-algorithms 
could  be  coded  and  tested  independently  within  the  main  framework  .  Figure  5.1 
shows  a  functional  flow  diagram  which  depicts  the  format  of  the  coding. 

Appendix  I  Is  a  listing  of  the  main  routine  and  the  subroutines  used  In  the 
final  algorithm.  Parts  of  the  code  use  system  unique  utility  programs  which 
were  developed  by  and  for  the  Magnavox  Image  Processing  Facility  and  are 
contained  In  the  system  library.  These  subroutines  are  principally  I/O  routines 
and  are  not  listed  In  the  appendix.  The  algorithm  code  was  written  In  a  fortran 
5  language  supported  by  the  Data  General  Eclipse  computer. 
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CONCLUSIONS  AND  RECOMMENDATIONS 


6.1  CONCLUSIONS 

A  good  detection  algorithm  must  have  a  high  sensitivity,  giving  a  high 
detection  rate,  but  It  must  also  have  a  high  rejectlvlty,  giving  a  low  false 
alarm  rate.  The  high  sensitivity  requires  the  ability  to  detect  low  wire 
signals  In  the  presence  of  high  noise  and  clutter  levels  which  corrupt  the 
signal.  This  means  that  It  must  In  some  way  average  out  the  noise  and  clutter 
fluctuations.  As  can  be  seen  from  the  video  tape  demonstration  accompanying 
this  report,  the  seml-llnear  detection  algorithm  developed  here  has  done  fairly 
well  In  this  area.  The  high  rejectlvlty  requires  the  ability  to  discriminate 
against  edges  and  points  as  well  as  against  line  segments  that  do  not  constitute 
extended  lines.  The  seml-llnear  detector  algorithm  does  quite  well  on 
discriminating  against  edges  and  points,  while  the  semi-local  line  discriminator 
algorithm  does  a  remarkable  Job  on  eliminating  Isolated  line  segir^  's  (line 
shaped  noise  or  clutter). 

Yet  It  appears  that  the  algorithm  Is  ttlil  no  substitute  for  the  human 
visual  system.  It  Is  evident  from  our  tests  that  the  detectivity  and  rejectlv¬ 
lty  of  the  human  eye/braln  system  when  interrogating  a  video  scene  Is  still 
superior  to  the  above  algorithm.  In  fact,  the  human  video  system  may  be  close 
to  the  limit  of  performance  that  can  be  expected  from  the  data,  particularly  In 
real  time.  Consider  that  the  human  eye  has  an  extraordinarily  wide  dynamic 
range  and  a  remarkable  sensitivity  In  this  range.  It  can  detect  differences  In 
average  brightness  of  as  little  as  IX  over  a  dynamic  range  of  nearly  1,000 
centered  about  the  eyes  adapting  brightness.  This,  taken  together  with  the 
eye/braln  ability  to  average  over  many  pixels  In  a  large  neighborhood  makes  for 
a  very  sensitive  detector.  Furthermore,  recent  evidence  exists  that  most  of  the 
cortical  cells  of  the  brain  respond  not  to  spots  but  to  specifically  oriented 


line  segments.  Add  to  this  the  brain’s  remarkable  ability  to  connect  local 
detections  In  one  area  of  a  scene  with  local  detections  In  another  area,  thereby 
giving  an  extraordinary  degree  of  rejection  of  detected  line  segments  that  do 
not  constitute  a  line  (e.g. ,  clutter),  and  one  has  an  extremely  efficient  line 
extracting  system. 

Nevertheless,  the  algorithm  developed  here  does  quite  a  creditable  Job 
and  It  may  be  superior  than  the  eye/braln  system  when  the  brain  Is  preoccupied 
with  a  variety  of  other  distracting  tasks  associated  with  flying,  target 
acquisition  and  fire  control,  etc.  The  algorithm,  as  It  stands,  may  be 
sufficient  when  combined  with  the  eye/braln  system  so  that  the  cable  enhanced 
video  gives  the  pilot  sufficient  warning  to  avoid  most  cables. 

On  the  other  hand,  there  Is  still  considerable  room  for  Improvement, 
and  there  are  two  areas  in  which  this  Improvement  could  be  accomplished.  These 
are: 

a.  Enlarge  the  local  neighborhood  sampling  to  Improve  the  noise 
filtering  effects  of  the  averaging.  This  area  certainly  would 
increase  the  computation  time,  but  should  also  result  In  a  more 
sensitive  detection  mechanism. 

b.  Develop  a  global  connection  algorithm  which  would  eliminate  the 
residual  clutter  noise  of  the  semi-local  line  discriminator. 

Again  this  would  Increase  computation  time  but  it  would  give  a 
degree  of  rejection  equivalent  to  the  eye.  It  may  be  possible 
here  to  Include  additional  features  for  discrimination  such  as  the 
constraints  of  hanging  cables.  This  could  lead  to  an  algorithm 
which  discriminates  against  lines  that  are  not  caused  by  hanging 
wires.  Remembering  that  every  increase  in  discrimination  ability 
allows  a  further  reduction  in  the  optimum  detection 


thresholds,  we  see  that  this  global  line  connection  algorithm 
would  lead  to  Increased  sensitivity  of  the  overall  system. 


As  noted,  these  Improvements  do  not  come  free.  They  will  require 
Increased  memory  capacity  and  Increased  computation  time.  The  current  version 
of  the  algorithm,  when  efficiently  coded  and  Implemented  In  hardware  should  be 
able  to  operate  In  real  time.  Only  a  thorough  Investigation  will  determine  the 
degree  to  which  these  algorithms  can  be  Improved  within  the  real-time  and  light 
weight  constraints  established  by  the  mission  requirements. 

6.2  RECOMMENDATIONS 

6.2.1  MAGIEC  feels  that  while  the  results  of  the  algorithms  developed  here 
are  very  encouraging  and  may  be  useful  as  they  stand,  it  Is  premature  to  enter  a 
hardware  Implementation  program  at  this  point.  MAGIEC  recommends  that  the 
hardware  development  effort  be  scheduled  to  commence  at  the  end  of  a  second 
phase  of  software  development.  We  feel  that  significant  Improvement  In 
performance  and  utility  may  be  gained  with  an  aggressive  software  Improvement 
program.  The  level  of  effort  expended  to  date  simply  was  not  sufficient  to 
develop  such  a  complex  algorithm  as  may  be  necessary  to  approach  the  performance 
of  the  human  visual  system.  This  software  Improvement  program  should  be  aimed 
at  those  areas  indicated  above,  namely: 

a.  Enlarging  the  neighborhood  sampling  of  the  detection  algorithm  to 
Increase  sensitivity. 

b.  Develop  a  global  connection  algorithm  to  discard  llne-llke 
clutter. 

We  believe  these  objectives  are  realisable,  but  they  are  sizeable  and 
will  require  a  substantial  level  of  effort  to  obtain  them. 

Once  a  global  connection  algorithm  has  been  developed.  It  may  be 
possible  In  the  future  to  develop  an  alarm  system  which  can  make  decisions, 
alert  the  pilot  to  potential  danger,  and  recommend  an  evasive  course. 
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IS  THE  MAIN  ROUTINE  USING  SAMSET 1 .  L.INCP2 .LINTES4 ,  THRE 
DIMENSION  NAi'iE(  5  ).  JJO(  9.0:255  >.MEIn:  5  >.  INAMEt  5  ) 
COMMON/RED/JDLK.LINE. J0<  0:255  ),JID( 4096  > 
COMMON/RIT/NDLK,LIN.KO( 0:255  )^1JQ( 4096  ) 

CDMMON/ANGLE/KA(  0:255  ) 

COMMON/U1NDOW/JV(  9,9  ).  JA(  3  >.  JB(  3  >.  JC(  3  ) 

COMMON/LINE/JVALC  B  ) 

JBLK=^K) 

NDLK^O 

LINE=0 

LIN=0 

DO  8  1=0.255 

8  K0(  I  )=KA(  I  >=0 

TYPE ’ENTER  READ  FILENAME" 

READ< 11.100  )NAME 
100  F0RMAT<5A2) 

OPEN  2.NAME.ATT="C' .REC=128 
OPEN  3."TEMP.DC* .ATT="C* .REC=256 
DO  9  1=0.3 

9  CALL  BRITE 

TYPE "ENTER  THRESHOLD  LINE  STRENGTH  * 

READ  FREEC 11 >ITHRES 
WRITE< 10.50  )ITHRES 
50  FORMAT<  "STRENGTH*" .15  ) 

C:ti^T:|;*READ  FIRST  EIGHT  ROWS5»5i|?#*it;;t: 

DO  1  J=2.9 
CALL  READl 
DO  1  JJ=0.255 
1  JJO(  J.  JJ  )=J0(  JJ  ) 

C*^:tc#DO  THE  WHOLE  PICTUREJ|!5KJI;:|:#* 

DO  2  K ROW =4. 251 

C-lJA^Ii^STORE  NINE  ROWS  FOR  THE  SAMPLING  MATRI X!|;:|c4:K'+ 

DO  3  K=1.8 
DO  3  J =0,255 
3  JJO(  K.  J  )=JJO(  K-ri.  J  ) 

CALL  READ! 

DO  4  J=0.255 
4  JJC<9.J>=J0<  J) 

C#*#*COMPUTF.  OPERATOR  FOR  THIS  ROW)»i5|;5l'>^l'.»SiS: 

DO  5  J=4.251 

CSi***COMPUTE  OPERATOR  FOR  THIS  PIXEL 
DO  6  K=l,9 
DO  6  L=1.9 

6  J9(K.L  )=JJ0<K.J+L-5  > 

DO  7  1=1.8 


CALL  SAMSETK  I  ) 

7  CALL  LIN0P2<  JVAL< I  ) ) 

CALL  LINTESKLINVAL.IVAL) 

IF<  IABS(  LINVAL  ).GT.6  )CALL  THRESK  ITHRES.LINVAL  .1 VAL  ) 
KA(  J  )=IVAL 
5  KO<J>=LINVAL 

CALL  DIRECT(KROW) 

2  CONTINUE 

DO  10  1=251.255 

10  CALL  BRITE 
CALL  RESET 

TYPE  ■<?>  <7>’ 

JBLK=0 

NBLK=0 

LINE=0 

LIN=0 

OPEN  2. 'TEMP.DC' .ATT='C' .REC=256 

OPEN  3,  NAME.ATT='C*,REC=128 

DO  11  1=0,255 

CALL  READ2 

CALL  BREAD 

DO  12  J=6.249 

IA.15,16 

14  K0<  J)=0 

15  GO  TO  12 

16  K0<  J  )=255 
12  CONTINUE 

CALL  QANTK* I.1,K0) 

11  CONTINUE 
CALL  RESET 
TYPE  ■<7> 

END 


V.,-1 


SUBROUTINE  SAMSETKN) 

C*!»:*9X9  sampling  CONFIGURATION. EIGHT  ORIENTATION' 
C0MM0N/WINri0U/J9<  9.9  ),  JA<  3  ).  JB<  3  )  .JC(  3  > 
JB(2>=J9(5,5) 

GO  TO  (10.20.30.40,50,60,70,80)  N 

10  JBC  1  >=J9(5,3  ) 

JB(3)=J9(5.7> 

JA<  1  >=J9<  3,3  ) 

JA(2)=J9(  3.5) 

JA(3)=J9<3.7) 

JC(1)=J9(7,3) 

JC<2)=J9<7,5) 

JC<3)=J9(  7.7) 

GO  TO  90 

20  JB(  1  )=J9(  6.3) 

JB<3  )=J9<4,7  ) 

JA(1  )=J9(  4,2) 

JA(  2  )=J9(  3,4  ) 

JA<3)=J9<2,6) 

JC(  1  )=J9(8,4  ) 

JC<2)=J9<7.6> 

JC(3  )»J9(  6,8  ) 

GO  TO  90 

30  JB(  1  )=J9<7,3  ) 

JB<3)=J9<3,7) 

JA(  1  )=J9(  5,1  ) 

JA(2)=J9(  3.3) 

JA(3  )=J9(  1,5  ) 

JC<  1  )=J9<  9,5  ) 

JC(2)=J9(7,7) 

JC(3)=J9<  5,9) 

GO  TO  90 

40  JB(  1  )=J9<  7,4) 

JB(3  >=J9<  3,6  ) 

JA(1  )=J9<  6,2) 

JA(2)=J9<4,3) 

JA(3)=J9(  2,4) 

JC<  1  )*J9<8,6) 

JC(2>=J9<6,7) 

JC(3)=J9(4,8) 

GO  TO  90 

50  JBC  1  )«J9(7,5  ) 

JB(3>«J9(  3,5) 

JA(1  )=J9(7,3) 

JA<2)»J9(5,3) 


JM  3  )=--J9(  3.3  ) 
JC(  1  )=J9<  7.7  ) 
JC<  2  )=J9(  5,7  ) 
JC(3  )=J9(  3.7  ) 
GO  TO  90 

60  JB(  1  )=J9(  7.6  ) 
JEK  3  >=J9(  3,4  ) 
JA(  1  )=J9(  8.4  ) 
J;.(  2  )=J9<  6.3  ) 
JA(3  )=J9(  4,2  > 
JC(  1  )=  J9(  6.8  ) 
JC(2  )=J9(  4.7  ) 
JC(  3  >=  J9(  2.6  ) 
GO  TO  90 

70  JB(  1  )=J9<  7.7  ) 
JB<3  )=J9(  3.3) 
JA(  J.  )=J9(  9,5  ) 
JA(2  )=J9(  7.3  ) 
JAv  3  )-=J9(  5.1  ) 
JC(  1  )=J9(  5.9  ) 
JC(  2  )=--J9(  3.7  ) 
JC<3  )=J9(  1.5  > 
GO  TO  90 

80  JB(  1  )=J9(  6.7  ) 
JB(  3  )=J9(  4.3  ) 
JA(  1  )~J9(  4.8  ) 
JA(  2  )='J9(  3,6  ) 
JA(  3  )=J9(  2.4  ) 
JC(  1  )=J9(  8.6  ) 
JC(2  )=J9<  7.4  ) 
JC(  3  )=J9(  6.2  ) 

90  RETURN 
END 


SUBROUTINE  LIN0P2(JVAL) 

C!Ic:»c*this  defines  the  line  characteristic  and/or  operator  to  be  used 

C0MM0N/WIND0U/J9(  9,9  ),  JA<  3  ),  JB(  3  >,  JC(  3  ) 

JASUM=JA(  1  )+JA(  2  )+JA<  3  ) 

JBSUM=  JB(  i  )+JB(  2  >+JB<  3  ) 

JCSUM=JC(  1  >+JC(  2  )+JC<  3  ) 

JUAL=0 

JBA=JBSUM-JASUM 

JBC=JBSUM-JCSUM 

IF(ISIGN< 1,JBA).EQ.ISIGN< l.JBC>>JVAL=JBA+JBC 
RETURN 


SUBROUTINE  LINTEBK LINUAL .1 UAL) 

C*5|;*THIS  COHPARES  ORIENTATIONS  AND  CHOOSES  THAT  WHICH  HAS  MAXIMUM  EVIDENCE 
COMMON/LINE/JVAL<  8  ) 

IVAL=0 
L.INVAL=0 
DO  10  1=1.8 
JV=IABS(  jOaL< I >) 

IF(  JV.GT.LINVAL  )1VAL=1 
10  IF< JV.GT.LINVAL)LINVAL=JV 
LINVAL=JVAL( IVAL  ) 

IVAL=ISIGN< IVAL.LINVAL  > 

RETURN 

END 


subrout; 

IF(  lABS 
IF<  lABSi 
RETURN 
END 


SUBROUTINE  DIRECKKROU) 

COMMON/RIT/NBLK, LIN. KO<  0:255  ).IJQ( 4096  ) 
COMMON/ANGLE/KA<  0: 255  ) 

COMMON/VALUE/KVAL<  5,0:255  ),KANG( 5.0:255  ) 
IF(KR0W.LE.7 )G0  TO  4 
DO  1  J=0,255 
DO  2  K=l,4 

KANG(K.  J  )=KANG(  K+1,J  ) 

2  KUAL(K,J)=KVAL(  K+l.J) 

KANG(5,J  )=KA(  J  ) 

1  KUAL<5,J)=K0(  J) 

DO  3  J=6.249 
K0(  J  )=KUAL<  3.J) 

IF<KO<  J).EQ.O)GO  TO  3 
CALL  CLUST(  J. JFLAG.nl ) 

IF(NL.LE.2)k6(  J)=0 
IF<  JFLAG.NE.2  )K0<  J  )=0 

3  CONTINUE 
DO  9  J=0.5 

9  K0(J)=0 
DO  10  J=250,255 

10  K0<  J>=0 
GO  TO  6 

4  DO  5  J=0.255 
KANG<  KROU-2 . J  >=K A<  J  ) 

KVAL(  KROW-2 ,  J  )=K0<  J  ) 

5  KA(  J  )=K0(  J  )=0 

6  IF(KR0W.GT.5)CALL  BRITE 
IF<KR0U.LT.251 )G0  TO  7 
DO  8  J=0,255 

8  K0<J)=0 
CALL  BRITE 

7  RETURN 
D 


SUBROUTINE  CLUSK J . JFLAG .NL) 
COMMON/VALUE/KVAL<  5,0:255  ),KANG(  5.0:255  > 
JFLAG=0 
NL=0 

K=IABS< KANG(3.J>) 

GO  TO  (10,  20.  30,  40,  50.  60,  70, ,80)  K 

10  KAP=KANG< 3. J+2  ) 

KAM=KANG( 3. J-2) 

GO  TO  90 

20  KAP=KANG< 2. J+2) 

KAM=KANG( 4. J-2  ) 

GO  TO  90 

30  KAP=KAN6( l,J+2 ) 

KAM=KANG<  5, J-2  ) 

GO  TO  90 

40  KAP=KANG<  1. J+1  ) 

KAM=KANG( 5. J-1  ) 

GO  TO  90 

50  KAP=KAN6<  l.J  ) 

KAM=KANG(5,J) 

GO  TO  90 

60  KAP=KANG(  1. J-1  ) 

KAM=KANG<5.J+1  ) 

GO  TO  90 

70  KAP=KAN6< 1, J-2  ) 

KAM=KANG( 5. J+2  ) 

GO  TO  90 

80  KAP=KANG<  2. J-2 ) 

KAM=KANG<  4. J+2  ) 

90  IF(  KAP.EQ.kANG( 3. J  ) )  JFLAG=JFLAG+1 
IF<  KAM.EQ.KANG<  3. J  )  >  JFLAG=JFLAG+1 
DO  1  K=-l.l 
DO  1  L=-l,l 

1  IF(  KANG<  3+K,  J+L).EQ.KANG<  3.  J  )  )NL=NL+1 
RETURN 
END 


END 


FILMED 


